A&A manuscript no. 

(will be inserted by hand later] 



Your thesaurus codes are: 
11.03.1, 11.06.2, 12.04.1 



ASTRONOMY 

AND 
ASTROPHYSICS 



0^ 



The fundamental plane of clusters of galaxies: 

a quest for understanding cluster dynamics and morphology 

Christoph Fritsch^ and Thomas Buchert^ '^ 

^ Mcix-Planck-Institut fur Extraterrestrische Physik, Giessenbachstrasse, D-85740 Garching, Germany 
^ Theoretische Physik, Ludwig-Maximilians-Universitat, Theresienstr. 37, D-80333 Miinchen, Germany 
^ Theory Division, CERN, CH-1211 Geneva 23, Switzerland 

Received 16 June 1998 / Accepted 3 December 1998 



(N 
> 
m 

o 

oo 

ON 
Oh 

6 ■ 



Abstract. We discuss implications of the fundamental 
plane parameters of clusters of galaxies derived from com- 
bined optical and X-ray data of a sample of 78 nearby 
clusters. In particular, we investigate the dependence of 
these parameters on the dynamical state of the cluster. We 
introduce a new concept of allocation of the fundamental 
plane of clusters derived from their intrinsic morpholog- 
ical properties, and put some theoretical implications of 
the existence of a fundamental plane into perspective. 

Key words: galaxies: clusters: general, fundamental pa- 
rameters; cosmology: dark matter 



1. The fundamental plane on galaxy scales 

The concept of the fundamental plane has its origin in 
some monovariant relations between global observables 
in stellar systems. The Tully-Fisher relation for spirals 
and the Faber-Jackson relation for ellipticals reveal that 
the orbital velocity increases strongly with the luminosity. 
However, the residual scatter suggests the introduction of 
an additional observable, the effective radius, which fi- 
nally leads to correlations in a three dimensional parame- 
ter space of observables, resulting in the concept of the 
fundamental plane (Dressier et al., 1987, Djorgovski & 
Davis 1987). Subsequently, this concept has been success- 
fully employed to assess the physical state of an elliptical 
galaxy. Observationally this plane is well established in 
the optical energy range (e.g. Guzman et al. 1993) and the 
infrared energy range (e.g. Pahre et al. 1995) separately, 
and is defined by a multivariate correlation between char- 
acteristic parameters of elliptical galaxies. Others like spi- 
ral galaxies have been shown to populate the same plane 
reaching down to the dwarfs by suitably adopting the pa- 
rameter space (e.g. Jablonka et al. 1996). 

Beyond the claim of existence of such a plane in some 
(observational, but not necessarily physical) parameter 
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space, it has proved to be a powerful concept to also inves- 
tigate the physical properties as well as the evolutionary 
state of the galaxy (e.g. Bahcall et al., 1995). Moreover, it 
can be employed as a distance indicator (e.g. Van Albada 
et. al., 1995). Thus, the fundamental plane has attained 
the status of a physical concept, i.e. a (conjectured) rela- 
tion among physical variables. 

2. Reaching out to the realm of cosmology 

It is a natural question to ask whether such a concept 
may be extended to larger spatial scales, but there are 
also strong implications of such an extrapolation which 
have to be seriously considered. 

2.1. The fundamental plane on cluster scales 

Although clusters of galaxies mainly consist of elliptical 
galaxies, this is not a straightforward argument to expect 
that such a plane can be found for structures on cluster 
scales. Schaeffer et al. (1993) have advanced that claim, 
however, based on a small set of 25 optically selected 
clusters. Recently, Adami et al. (1998) have consolidated 
this concept with a sample of 29 clusters within the ESO 
nearby Abell cluster survey. In the present note we advo- 
cate this claim. However, we would like to stress that the 
luminous matter in clusters of galaxies is mainly built from 
the X-ray emitting intergalactic gas (10% - 30%) and not 
from galaxies (1%-10%); the allocation of a fundamental 
plane should therefore be decided on the basis of X-ray 
data rather than optical data alone. Their combination 
carries additional information about the coupling between 
galaxies and intergalactic gas, which are both supposed to 
be settled down in the gravitational potential built from 
the hitherto unspecified 'dark matter' (70%-90%). Based 
on the work by Fritsch (1997) and Fritsch & Bohringer 
(in prep.), which supports the existence of a fundamental 
plane in the characteristic optical and X-ray properties of 
galaxy clusters, we are going to investigate implications of 
morphometric parameters that quantify cluster substruc- 
ture. 
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There is also a dynamical reasoning for the existence of 
a fundamental plane which assumes the scalar virial theo- 
rem, which was derived for isolated systems, to hold. This 
idealization provides a relation between the relevant phys- 
ical parameters, here: the gravitational mass, the extent of 
the structure under consideration (measured in terms of 
half-light radius, see below) and the stabilizing dynamical 
pressure due to velocity dispersion. We shall also consider 
this relation below to derive values for the mass-to-light 
ratio. Although simple-minded, the use of this relation 
stems from the common implication that an ensemble of 
"virialized" entities should define the fundamental plane. 
Can this be so simple? 

2.2. Theoretical implications 

While galaxies are easily identified as individual entities, 
although their dark halos may not fully follow that iden- 
tification, the cluster as a structural unit is more difficult 
to assess. People tend to consider clusters of galaxies as 
the largest gravitationally bound systems in the Universe, 
and the wordings "decoupled from the universal expan- 
sion" and "relaxed, or virialized system" are applied to 
suggestively establish the possibility of isolation from the 
large-scale structure environment. It is clear that this can 
only be true for an idealized cluster; real ones are nei- 
ther "relaxed" , nor "decoupled" from the expansion and 
their environment. To write down the simplest possible 
formula for the relation between physical variables (the 
scalar virial theorem), and to use this relation as a stan- 
dard of reference of a "relaxed" cluster is, at least, coura- 
geous. 

What is a cluster? 

In general we may consider any overdense patch of matter 
and require some stationarity condition to hold for it on 
a,verage. This defines the term "dccoiipling from the ex- 
pansion" . The tensor virial theorem as defined by Chan- 
drasekhar & Lee (1968) cannot be applied per se to a 
non isolated system. Instead a generalized cosmic virial 
theorem should state the relation between spatially av- 
eraged quantities for some spatial domain embedded into 
the cosmological model. This, still, requires the treatment 
of boundary terms that arise by averaging over an over- 
dense portion of matter at the boundaries of the averag- 
ing domain. Within the cluster all dynamical variables in 
general fluctuate, and in turn the strength of the fluctua- 
tions themselves has influence on the average properties. 
The latter is known as the "backreaction effect" (see, e.g., 
Buchert & Ehlers 1997 and ref. therein). The picture to 
be emphasized here is that a cluster which is embedded 
into the large-scale structure environment should be sub- 
jected to a stationarity condition on average taking the in- 
ternal dynamics and cluster boundary terms into account 
(compare here the so-called C-correction that has been 
recently applied by Girardi et al. 1998). 



The m,erging problem,: 

Dynamically, since we face the problem of merging, we 
may pick that specific clump of matter belonging to a sta- 
tionary system and trace this matter back into the past by 
conserving its mass (the Lagrangian point of view) . We so 
avoid that there is matter in- or outflow accross the clus- 
ters' boundary. The evolution of morphological properties 
of this clump of matter may then be considered, map- 
ping the dynamical state as a function of time in terms 
of morphological parameters. It is here, where the idea of 
using morphological descriptors may provide a landmark 
of how to diagnose the parameter space in which we want 
to allocate a fundamental plane. 

The morphology-cosmology connection: 
The study of the evolution of cluster morphology is a lively 
debated subject in connection with cosmological simula- 
tions. Especially the relation to the background cosmol- 
ogy is the focus of interest in this field (Evrard et al. 1993, 
Mohr et al. 1995, Crone et al. 1996). Still, we consider it 
important to understand the notion of a "relaxed" clus- 
ter from first principles. If we talk about the existence 
of a fundamental pkme, then wc imply this only for some 
asymptotic dynamical state of "relaxed" clusters for which 
the relation implied by the fundamental plane is exactly 
satisfied. This relation must be sought theoretically. We, 
here, see an intimate link between the dynamics of the 
cluster and a conjectured attractor in the space of the 
characterizing averaged variables. 

We suggest to quantify the deviation from these yet 
unknown defining properties of "relaxed" clusters by in- 
strinsic morphological properties such as their amount of 
substructure. We explain this idea now for a sample of 
clusters based on optical and X-ray data. 

3. The optically selected and X— ray based cluster 
data 

3.1. The sample 

Using the COSMOS-/ APM- and the ROSAT X ray data 
robust optical and X-ray parameters for a subset of 78 
clusters of galaxies from surface brightness profiles were 
derived in (Fritsch 1997; for details concerning the data 
see Fritsch & Bohringer, in prep.). To avoid the influence 
of evolutionary eff'ects a homogeneous sample of clusters 
of galaxies in a redshift range of 0.02 < z < 0.05 was 
studied. A family of parameters was so derived that, at 
first phenomenologically, characterize the physical state of 
the clusters. The most important independent parameters 
are the luminosities and the half-light radii Rq (optical) 
and Rx (X-ray). 

Fitting ellipses to each of the projected distributions 
gives us the optical/ X-ray centers of the clusters. The cal- 
culation of the background is based on the local galaxy/ 
photon distribution outside the cluster. Background cor- 
rected, differential surface brightness profiles allow us to 
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measure the radii that contain the total Ught of the clus- 
ters. Finally, the half light radii were calculated from in- 
tegrated surface brightness profiles. They define the radii 
within which half the light of the cluster is emitted. The 
total light Lq of the clusters emitted from the galaxies 
within the optical blue band is given by all the background 
corrected galaxy magnitudes within a circle that contains 
all the light. To calculate the light below the given mag- 
nitude limit we apply the Schechter luminosity function 
with = - 21.8 and a = -1.25. 

The total X-ray luminosity was calculated iteratively 
by using the Raymond-Smith code (Raymond & Smith 
1977) for a completely ionized plasma with a typical 
metallicity of half the solar one, and the empirical corre- 
lation between the X-ray luminosity and the temperature 
kT oc L"-^^'* given by White (1996). The latter correlation 
was used because ROSAT does not provide a reasonable 
determination of the temperature and the literature offers 
too few X-ray temperatures for a correlation analysis as 
far as our clusters are concerned. Superpositions of other 
sources were detected via cross correlation with clusters 
stored in the NED-database (NASA/IPAC Extragalactic 
Data Base) and cut out, subsequently. 

For all the uncertainties of the derived parameters the 
errors due to spatial binning and the errors of the source 
and background count statistics were included. 

The complete list of the 79 clusters is stored in Table ||. 

3. 2. Morphological method of allocation 

That clusters of galaxies are not arbitrarily distributed in 
the three-dimensional parameter space {Lq, L^, Rq} is not 
surprising, but that they lie within a fairly well-defined 
plane is enough reason to introduce the concept of the 
fundamental plane of clusters. To find out the physical 
meaning of this phenomenological "fundamental plane" 
is another issue. As outlined above we may approach the 
problem from first principles. Here, we would like to sketch 
a procedure which already lays down a fairly unique way 
to establish a diagnostic criterion of allocating the fun- 
damental plane. We stress, however, that in practice the 
sample of clusters has to be larger than the one we are 
going to study in order to get statistically significant re- 
sults. 

Already from visual inspection we appreciate that the 
clusters can be classified in terms of some intrinsic struc- 
tural property: the clusters should pass the test whether 
they are useful to define the fundamental plane in a phys- 
ical sense. We have to find a measure of the dynamical 
state of the clusters of galaxies. Considering the process 
of structure formation and the relaxation process for the 
galaxies and the g dynamical state may be repre- 
sented by some morphological measure which character- 
izes the amount of substructure in the clusters. We adopt 
the working hypothesis that "relaxed" clusters of galax- 
ies show a small amount of substructure and "unrelaxed" 



ones, which are still in the process of merging, show strong 
substructure. Useful substructure measures have already 
been proposed and employed in the literature (e.g.. Crone 
et al. 1996 and ref. therein). We, here, base our measure 
on radial variations of morphological parameters derived 
from fitting ellipses to the projected distribution of the 
galaxies and X-ray photons which includes the center-of- 
mass shift, the ellipticity and the position angle of clus- 
ter contours (Fritsch 1997). In this line, robust structure 
functions based on vector-valued Minkowski functionals 
have been proposed recently (Beisbart & Buchert 1998) 
and are currently tested on simulated clusters. This new 
method will also help to overcome possible biases in the 
presently used morphological method which does not dis- 
tinguish substructure from "twisted" isocontours. 

3.3. Laying down the fundamental plane of clusters 

Taking these substructure measures we can divide the 
whole sample of clusters due to the amount of their sub- 
structure into two classes with the same number of mem- 
bers. Now we consider the polynomial Pep approximat- 
ing the data for Lq: PEp{Ra,L^) = E^j^^ l0-2i (Fig. 
1). It results from a two-dimensional x^-fit according to 
the Levenberg-Marquart algorithm to all the cluster data 
in the three-dimensional parameter space, consisting of 
the optical luminosity Lq, the X-ray luminosity and 
the optical half-light-radii i?o (we call this the empirical 
plane of clusters) . The orthogonal scatter of that plane is 
given by 24 %. Applying a correlation analysis for both 
classes separately between the optical luminosity Lq and 
the function Pep{Ro, L^) we determine the probabilities 
for the null hypothesis that there is no correlation between 
Pep{Ro, Lx) and Lq; we obtain P(rp^p:io) ^ 10^^ for the 
clusters with less substructure, and P(Tp^p;i) ^ 10^^ for 
the clusters with much substructure^. Therefore, the for- 
mer class may serve as that sample which phenomenologi- 
cally comprises the "relaxed" clusters (defining the funda- 
mental plane) with a reduced orthogonal scatter of 21 %. 
Additionally we find a strong correlation between the dis- 
tance of the clusters to the so-defined fundamental plane 
and the corresponding substructure measure (Figs. 1,2). 
This correlation supports the hypothesis that the location 
of the fundamental plane is related to the feature of 'less 
substructure'. 

Whether this morphological property uniquely relates to 
the virial condition, and whether the location of the funda- 
mental plane is refiected correctly by the standard use of 
the virial relations among the parameters is not clear at all 
(see also Adami et al. 1998). However, in order to infer in- 
formation on the mass-to-light ratio, we have to apply the 
virial condition in its usual, albeit naive form: to estimate 
the masses of the clusters of galaxies we use the empirical 



Tp^p:La is Kendall's t quantifying the correlation between 
Pep{Ro, Lx) and Lo in a non-parametric way. 
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relations between velocity dispersion and temperature in 
consistency with the thermodynamical equilibrium condi- 
tion (cr^ oc T) , and between temperature and X-ray lumi- 
nosity (T oc L°-^^^) given by White (1996). Additionally, 
we have to assume a density profile to find the relation be- 
tween the optical half-light radius and the gravitational 
radius Rg . For simplicity we model the distribution of the 
whole matter with a King profile; note, however, that the 
cluster mass may be over-/underestimated by a factor of 
a few, if the true density profile is steeper/shallower in 
the core of the cluster (see: Sadat 1997). Then we start 
from the observed correlations between the parameters 
and our morphology-based allocation of the plane and 
transform these parameters into the physical parameter 
space M,a,Ro (Fig. 3). 

3.4- Mass-to-light ratios for clusters of different 
morphology 

Given the observed relations and the standard virial con- 
dition, = ^T^, the empirical plane based on the total 
sample can be represented by the characteristic mass to 
optical light relation: 

f M y f yf^Q\ 

with coefficients listed in Table |l| under 'EP'. 

From Table ^ we infer that there is a clear trend of 
discrimination between clusters with much substructure 
(index 'SP') and clusters which belong to the class of more 
"relaxed" clusters (i.e. with less substructure in the opti- 
cal and X-ray energy range) (index 'FP'). However, the 
respective data sets still overlap within the errors which 
is a consequence of the very small sample of clusters on 
which we base our analysis. Even if we don't admit any sig- 
nificant difference between the empirical and fundamental 
planes, there is a trend that the characteristic mass-to- 
light ratio for the "relaxed" clusters depends only slightly 
on the mass and almost does not depend on X-ray lumi- 
nosity; it suggests a constant value for M/L which is a 
striking result of our morphology-based allocation of the 
fundamental plane. The latter would also support the as- 
sumption that the light distribution of the galaxies follows 
the mass distribution in "relaxed" clusters, an assumption 
that lies on the basis of most mass estimates (see, e.g., 
Girardi et al. 1998). Furthermore, in comparison with the 
mass-to-light ratio for elliptical galaxies (Bender et al. 
1992), our data imply that the mass-to-light ratio of clus- 
ters reveals about 30 times more "non-optical luminous 
mass" than the elliptical galaxies which gives a hint to a 
certain fraction of the X-ray mass and especially to the 
amount of the underlying unknown dark matter (Fig. 4). 
It is interesting that our analysis suggests a tendency to- 
wards lower values for M/L in clusters as we approach the 
fundamental plane. 



Table 1. Coefficients for the mass-to-light ratios for the 
subsamples 'EP' (empirical plane, all clusters), 'FP' {fun- 
damental plane), 'SP' (clusters with much substructure). 



index 


a 


b 


c 


EP 


448.70 ± 148.21 


0.16 ±0.17 


0.09 ±0.12 


SP 


613.16 ±190.17 


0.30 ±0.14 


0.08 ±0.15 


FP 


318.32 ± 107.13 


0.12 ±0.18 


0.10 ±0.13 



4. Conclusions 

The data show that the nearby (0.02 < z < 0.05) clusters 
of galaxies lie preferentially within a plane in a three- 
dimensional parameter space built from the optical lumi- 
nosity, the X-ray luminosity and the optical half-light ra- 
dius. On the basis of a morphological criterion for cluster 
substructure and the hypothesis that "relaxed" clusters 
reveal less substructure, we identified a fundamental plane 
by a 2-dimensional fit to a (by a factor of 1/2) reduced 
sample of clusters with less substructure. The distances of 
the other clusters to this fundamental plane show strong 
correlations with our measure of substructure. The pro- 
posed method for allocating the fundamental plane should 
be viewed in parallel with a physical condition among spa- 
tially averaged dynamical variables, which still has to be 
sought. Adopting the generally held view of the standard 
virial relations (based on the trace of the tensor virial the- 
orem for isolated systems), we derived the M/L values for 
each class of clusters implying a tendency towards lower 
values for more "relaxed" clusters. Our method suggests 
a weak dependence on mass and luminosity for clusters 
populating the fundamental plane resulting in an almost 
constant value for M/L of typically 300. The (orthogo- 
nal) scatter around the empirical plane of 24% for the 
sample of all 78 clusters is reduced to a scatter around 
the diagnostically selected fundamental plane of 21 % for 
the sample of 29 clusters with less substructure. Both are 
notably amplified when the L^ dependence is ignored. 
For future work, this scatter can be analyzed in more de- 
tail, if one takes into account different heating mechanisms 
for the intergalactic gas, different profiles of the intergalac- 
tic gas and the dark matter, and different populations of 
galaxies within single clusters. Furthermore, the scatter 
around the planes may be used as an indicator of the evo- 
lution of clusters, if one applies the concept of the funda- 
mental plane to samples of clusters belonging to different 
redshift ranges. 
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Fig. 1. The empirical plane for clusters of galaxies fitted to the whole sample in face-on projection (shown in the 
log(Lo/L0) — log(_L2;/10''"'erg/s) — log(i?o)-space). The grid-plane results from a 2-dimensional fit to the data. The 
black symbols mark clusters of galaxies with less substructure in the optical and X-ray energy bands (members of the 
fundamental plane). The white symbols mark the complementary sample of clusters with more substructure. 
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Fig. 2. The empirical plane (Lq = F{^'^ -^2'^^) ^oi clusters of galaxies fitted to the whole sample in edge-on projection 
(marked by the line with steeper slope). Again, the black dots mark the clusters with less substructure, whereas open 
circles mark clusters with more substructure in both the optical and X-ray energy bands. The line with shallower 
slope represents the fit to the more "relaxed" clusters defining the fundamental plane 
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Fig. 3. The velocity dispersion versus virial mass. The velocity dispersion for the clusters is derived from the measured 
X-ray luminosities using the empirical relations between the X-ray luminosity and the temperature kT oc given 
by White (1996). The circles denote our clusters (including a white dot, if belonging to the fundamental plane). The 
squares denote the clusters with effective radii taken from Cappi & Maurogordato (priv. comm.) and velocity dispersion 
taken from Struble & Rood (1991). The line with steeper slope represents the fit to all clusters corresponding to the 
empirical plane. The other line represents the fit to the more "relaxed" clusters defining the fundamental plane. 
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Fig. 4. The empirical plane fitted to the whole sample in the log(Lo/i0) ~ log(Ai"/Ai^0)-projection for structures on 
different scales. Black circles mark clusters of galaxies and white squares mark elliptical galaxies taken from Bender 
et al. (1992). The upper line corresponds to objects, which would consist only of optically luminous matter. 
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Table 2. List of the clusters with the luminosities Lq (optical) and (X-ray), and the optical half-light radii Rq. 



cluster 


Lo 




Rq 






fl0*^erg/sl 

L ^ o/ "J 


[Mpcl 


a0076 


143.074 


1.803 


0.833 


a0119 


451.292 


11.061 


1.559 


a0147 


157.714 


1.462 


1.482 


a0160 


220.267 


0.949 


1.240 


a0168 


210.932 


2.665 


1.286 


a0189 


99.819 


1.340 


2.098 


a0195 


134.423 


0.307 


0.413 


a0260 


290.824 


1.294 


1.457 


a0261 


12.599 


0.067 


0.247 


a0295 


43.532 


0.653 


0.556 


a0376 


233.546 


4.204 


0.939 


a0400 


25.010 


1.233 


1.394 


a0407 


199.971 


2.302 


1.607 


a0533 


124.108 


0.712 


1.030 


a0576 


172.687 


4.679 


0.932 


a0779 


85.412 


0.345 


1.188 


a0999 


82.340 


0.136 


1.290 


allOO 


121.335 


0.688 


0.940 


all39 


123.110 


0.706 


1.138 


all42 


72.391 


0.720 


1.170 


all77 


74.972 


0.513 


1.645 


all85 


153.499 


0.938 


1.193 


al213 


151.170 


0.925 


1.088 


al228 


241.229 


0.266 


1.165 


al314 


201.711 


0.990 


1.063 


al367 


157.068 


3.069 


0.705 


al644 


296.538 


9.702 


0.720 


al656 


487.000 


8.000 


1.409 


al736 


322.886 


8.790 


1.150 


al983 


410.034 


1.237 


1.713 


a2052 


122.214 


8.781 


0.600 


a2063 


439.425 


6.862 


2.711 


a2107 


120.115 


1.755 


0.379 


a2147 


306.000 


10.023 
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